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摘 要 :针对 建筑 结构 基于 Clough-Penzien(C-P) 谱 下 的 结构 随机 地 震动 响应 既 有 方法 复杂 的 现状 ， 

提出 了 一 种 简明 解析 法 。C-P 谱 是 将 场地 的 地 震动 看 作 基 于 白 骂 声 激励 的 两 个 关联 的 滤波 方程 , 较 
准确 地 模拟 场地 的 平稳 地 震动 特征 。 首 先 ,提出 将 C-P 谱 滤波 方程 组 与 建筑 结构 的 运动 方程 联 立 ， 
实现 了 不 易 求 解 的 非 白 嗓 声 激励 精确 转化 为 简明 的 白 噪 声 激励 的 效果 ;其 次 ,基于 复 模 态 正 交 法 对 
重 构 后 的 结构 运动 方程 组 进行 解 耦 ,并 基于 随机 振动 理论 推导 出 结构 响应 的 功率 谱 密度 函数 和 0 
阶 、1 阶 和 2 阶 谱 矩 的 简明 解析 式 ; 最 后 ,分 析 了 基于 首次 超越 破坏 准 和 最 弱 失 效 模式 下 的 结构 体系 
动力 可 靠 度 。 算 例 验 证 了 所 提 方 法 的 准确 性 和 简洁 性 ;此 外 ,所 提 方 法 也 可 推广 到 其 他 复杂 平稳 非 
白 骂 声 激励 的 地 震动 响应 分 析 。 
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Abstract: Aiming at the inefficiency of structural response and dynamic reliability analysis under Clough- 
Penzien (C-P) random earthquake excitation ,a simple and fast method for solving the seismic response e- 
quation of building structure by using C-P filtering equations and seismic response equation of building 
structure is proposed. C-P spectrum regards the ground motion of site as two filtering differential equa- 
tions, which can accurately simulate the stationary ground motion characteristics of site. First,the complex 
ground motion excitation is transformed into the ground motion excitation based on white noise by combi- 
ning filtering equations of C-P spectrum with structural dynamic equations. Second, the complex mode 
method is used to decouple the ground motion equations ,and the concise analytical expressions of covari- 


ance, power spectral density function and 0-2 order spectral moments of structural response are obtained 
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based on the random vibration theory. Last, the system dynamic reliability calculation of the structure 


based on the first passage failure criterion is analyzed. An example shows the accuracy and simplicity of 


the proposed method. In addition ,this method can also be applied to the analysis of ground motion re- 


sponse subjected to other complex stationary non-white noise rand excitation. 


Key words:Clough-Penzien (C-P) spectrum ;spectral moment;system dynamic reliability ; concise analyt- 


ical solution; complex mode method 


地 震 时 常 发 生 ,每 次 强 震 都 会 导致 震 区 一 定 范 
围 的 各 种 人 造 结构 诸如 房屋 建筑 ,交通 设施 .水利 设 
施 等 的 损毁 并 造成 不 同 程度 的 人 员 伤亡 。 因 此 ,对 
于 地 震动 的 研究 由 来 已 久 , 自 20 世纪 40 年 代 开始 ， 
随 着 地 震动 观测 设备 的 改进 及 观测 数据 的 增加 ,发 
现 地 震 具 有 显著 的 随机 性 ,各 种 随机 地 震动 模型 被 
提出 。 文 献 [1] 首先 提出 白 噪 声 激励 模型 ,该 模型 
认 态 功率 谱 在 全 频率 范围 内 保持 不 变 ,这 与 真实 的 
地 震 过 程 存 在 偏差 ,但 开启 了 随机 模型 描述 地 震动 
的 铝 河 。 文 献 [2-3] 根 据 大 量 的 地 震动 时 程 记 录 , 提 
出 忆 过 滤 白 噪声 模型 , 即 著名 的 Kanai-Tajimi 随机 


算 精 度 受 数 值 积 分 参数 取 值 的 影响 。 

时 域 法 ,以 激励 的 协 方差 为 基础 ,以 结构 的 结构 
响应 方差 为 研究 对 象 , 实 模 态 和 复 模 态 是 其 两 种 主 
要 方法 "2 。 实 模 态 方法 仅 能 解决 经 典 结构 的 振 
动 问题 ,而 复 模 态 方法 是 一 种 解决 定常 高 阶 微分 方 
程 的 数学 方法 ,因而 可 广泛 应 用 于 各 类 线性 结构 的 
振动 问题 ”““ 。 文献 [7] 研 究 了 Clough-Penzien 谱 
下 隔 振 结构 的 随机 响应 ,所 获得 的 响应 表达 式 比较 
复杂 , 且 仪 给 出 了 结构 的 方差 ,没有 给 出 工 阶 谱 矩 的 
分 析 。 

基于 首次 超越 破坏 准则 结构 地 震动 可 靠 度 分 析 


谱 * 该 模型 认为 场地 土 下 覆盖 的 基 岩 的 地 震动 过 程 
汶 自 噪声 过 程 ,而 覆盖 士 层 的 运动 看 作 一 考虑 场地 


是 结构 安全 评估 的 重要 手段 5 , 目前 工程 上 采用 
以 结构 层 间 变形 的 失效 模式 为 弹性 结构 的 失效 模 


上 十 特 征 的 单 自由 度 过 滤器 ,能 解决 白 噪 声 激励 模型 
功 雍 谱 变 化 的 问题 ,但 存在 过 分 夺 大 低频 地 震动 能 
量 的 缺点 。 随 着 地 震动 记录 的 不 断 增加 ,各 国学 者 
均 提 出 了 新 的 地 震动 模型 ,如 美国 学 者 Clough 和 
Fenzien 提出 的 Clough-Penzien 谱 呈 中 国学 者 胡 奸 
贤 敬 提出 的 胡 替 贤 谱 ”、 欧 进 萍 等 提出 的 欧 进 薄 
谱 过 等。 其 中 , Clough-Penzien( C-P) 谱 利用 2 个 滤 
小 圳 来 描述 场地 土 对 地 面 运动 的 影响 , 可 改善 Ka- 
nafTajimi 谱 不 足 , 在 工程 上 有 着 广泛 的 应 用 "? 。 
随机 地 震动 下 建筑 结构 的 响应 分 析 法 主要 有 时 
域 法 和 频 域 法 “， 。 频 域 法 应 用 时 ,结构 响应 功率 
谱 表示 为 地 震 激励 功率 谱 与 结构 的 频 响 函 数 模 值 平 
方 的 代数 积 ,无 需 解 微分 方程 ,而 受到 广泛 的 应 用 。 
最 具有 代表 性 的 实用 方法 为 虚拟 激励 法 "" 和 传递 
函数 法 ,其 中 ,林家 浩 提出 的 虚拟 激励 法 已 广泛 
应 用 于 各 种 结构 的 振动 分 析 研 究 ，”。 李 春 祥 
等 "利用 传递 函数 法 研究 了 基于 修正 C-P 谱 的 
MTMD 结构 的 动力 特征 。 李 创 第 等 ' 利用 传递 函 
数 法 研究 了 基于 巴 斯 金 谱 的 设置 支撑 的 广义 Max- 
well 阻尼 器 系统 的 地 震 响 应 分 析 。 然 而 工程 中 需要 
计算 结构 随机 响应 的 方差 和 谐 矩 等 参数 ,以 上 研究 
的 参数 频 域 法 都 需要 在 频 域内 进行 数值 积分 , 对 于 
大 型 复杂 结构 的 随机 响应 分 析 , 计 算 效 率 比较 低 , 计 


式 ,采用 层 间 弹性 位 移 角 作为 指标 。 对 于 多 自 
由 度 体系 的 动力 可 靠 度 分 析 , 首 要 前 提 是 确定 其 失 
效 模式 , 欧 进 萍 等 ”提出 了 “最 弱 失 效 模式 ” ,得 到 
了 广泛 地 应 用 ” 。 刘 强 等 ” 利用 虚拟 激励 法 研 
究 了 小 变 结 构 体 系 基于 Poisson 分 布 的 动力 可 靠 度 ， 
该 研究 所 得 的 位 移 方差 和 速度 方差 均 为 数值 解 ,其 
分 析 效 率 及 精度 受 数 值 积分 间隔 影响 较 大 。 

本 文 将 Clough-Penzien 谱 的 2 个 滤波 方程 与 结 
构 的 运动 方程 联 立 ,将 复杂 的 地 震动 激励 准确 地 转 
化 成 简明 的 白 噪 声 激励 ;基于 复 模 态 正 交 法 和 随机 
振动 理论 ,推导 结构 随机 地 震动 响应 的 协 方差 谱 和 矩 
等 响应 量 的 简明 解析 式 ; 最 后 ,基于 首次 超越 破坏 准 
则 和 “最 弱 失 效 模式 "分析 结构 的 体系 动力 可 靠 度 。 


1 结构 地 震动 方程 基于 C-P 谱 的 白 
噪声 转换 


一 nn 层 建筑 结构 ,其 质量 .阻尼 、 刚 度 矩 阵 分 别 

为 MC、K, 地 震动 下 的 运动 方程 为 
MX+Cx+Kx= -MIX, (1) 
式 中 :x* xx 分 别 为 建筑 结构 相对 于 地 面 的 加 速度 、 
速度 和 位 移 ;T = [1 1] ,上 标 “T” 表 示 向 量 转 
置 运算 ;x, 为 场地 地 面 运 动 的 绝对 加 速度 ,本 研究 
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采用 C-P 谱 的 随机 地 震动 模型 , 该 模型 的 滤波 方 
程 为 “ 


tg 26rorxs + Wg = 一 %1 (2a) 
Xj=Us tunp (2b) 
Ls +2 WU tO = — ip (2c) 


式 中 :w, .é, 分 别 为 基 岩 覆盖 层 上 场地 土 的 阻尼 比 和 
车 越 频 率 ;wj.é 分 别 为 地 面 处 场地 土 的 阻尼 比 和 卓 
越 频率 ; ii, 为 场地 处 基 岩 地 震动 加 速度 ,为 白 噪 声 
随机 过 程 ,其 协 方差 函数 为 
Ci (7) =2mS08(7) (3) 
式 中 : 5。 为 基 岩 地 震动 强度 常数 ; 8 (r) 为 Dirac 
函数 。 
工程 上 常 采用 C-P 谱 的 功率 谱 密度 函数 ”0 , 即 


ka) = ww ow +4é ww 
3 0) ro (0 0) toe 
(4) 


中 由 式 (4) 可 知 ,区 的 功率 谱 密度 函数 表达 式 极 
其 受 杂 ,在 工程 应 用 时 无 法 获得 简明 的 解析 解 。 为 
些 f 杀 文 提出 了 基于 C-P 谱 激 励 模 型 下 结构 的 地 震 
动 痊 析 可 采用 其 过 滤器 方程 与 建筑 结构 的 地 震动 运 
动 刘 程 联合 建立 分 析 的 方法 。 

之 联 立 式 (1) 及 式 (2) , 则 建筑 结构 的 运动 方程 的 
矩阵 形式 为 


MX+CX+RIY= -Qi (5) 
X= [us Xs x] (6a) 
1 0 0 
M=I1 1 0 
0 MI Mj,,,, 
w 0 0 
Ee (6b) 
0 0 K n+2,n+2 
2é 0, 0 0 1 
| 0 2é 0 = 
C n+2,n+2 0 | 
(6c) 


由 式 (5) 可 知 ,建筑 结构 基于 C-P 谱 的 地 震动 
转化 为 白 噪 声 激励 的 结构 地 震动 ,为 本 文 所 获得 C- 
P 谱 下 地 震动 响应 的 简明 解 英 定 基础 。 

2 建筑 结构 地 震动 响应 的 统一 解 


引入 状态 变量 


和 


y=[X XX] (7) 
式 (5 ) 改 写 为 
My +Ky=run(t) (8) 
式 中 
-~ IM 01, 10 -天 0 
的 | 人 =| C = -| 


依据 复 模 态 理论 "…”"! ,方程 式 (8) 存 在 使 其 解 耦 
对 角 阵 己 及 左右 特性 向 量 UV, 即 
V' KU 


(9) 


y=Uz (10) 
式 中 z 为 广义 复 模 态 变量 。 
利用 式 (10 ) 并 结合 复 模 态 正 交 法 的 特性 , 式 
(8) 简 写 为 


z+Pz=m un (11) 
式 中 
Vr 
了 vv 
由 于 P 为 对 角 阵 , 可 获得 式 (11) 的 分 量 式 
Z; +p;z; = Nunp (12) 


式 中 zi\p; 分 别 为 向 量 z,m,p 的 分 量 。 

联 立 式 (6a) 式 (7) 式 (10) 及 式 (12) ,结构 层 
的 位 移 x; 和 速度 x;, 层 间 位 移 Ax; 及 层 间 位 移 变 化 
率 Ax;, 可 表示 为 


2(n+2) t 


wv =u = 子 Rial it -7)dr (13a) 
j= 0 


(n+2) 


2 1 
a pir 
Yi Ujirn+42 >» havin LRC 7) dr 


i=1 0 


(13b) 
2(n+2) v 
Ax; = 之 和 -7T)dr,] > 1 
(13c ) 
2(n+2) » 
A = 也 Vy | al -7)drj > 1 
= 0 
(13d) 
AN = Ws A = (13e) 
结构 响应 的 强度 系数 为 
Aji = wim; (14) 


式 中 心 ,as = (52 -wi) i 
结构 的 各 类 响应 可 统一 表示 为 


2(n+2) 上 2(n+2) 


Ny Dx fe ialt -7) dr = EX) 


(15) 
式 中 为 结构 响应 的 强度 系数 。 
结构 的 响应 分 量 为 
XD) =x fe ilt-r)dr (16) 


3 随机 响应 功率 谱 和 谐 矩 分 析 


3.1 随机 响应 的 协 方差 分 析 

依据 随机 振动 理论 , 并 通过 式 (15 ) ,结构 响应 
X 的 协 方差 表示 为 
有 2(n+2)2(n+2) 
人 cx) = 2 SELXAOX(t + 7)] (17) 
由 式 (16) ,随机 响应 分 量 的 协 方差 表示 为 
X(t + 7)] 
TA, 下 | eewElial: —u) ia(t +7T -v0)]dudv 


Qt [ | eeeeC (uw +7 -v0) dudy (18) 
E 0 J0 a 


PL X(t) X(t + 7)] 


A i,x | | ered(n a 
克利 有 Dirae 的 函数 性 质 , 式 (19) 简 化 为 
EC 4 7)] = 2 xx, | me 
(20) 
对 式 (20) 运 用 积分 
ELX(DOXC+7)] 275 XX (21) 


由 式 (17) 及 式 (21) ,结构 基于 C-P 谱 的 地 震动 
响应 的 协 方差 为 


2(n+2) 
Cx(7) = 2750 > Gue™” (22) 
el 
2(n+2) Xx 
人 二 AAA (23) 
外 名 py 十 万 


观察 式 (22 ) ,建筑 结构 基于 C-P 谱 的 随机 地 震 
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2(n+2) 
ox = VCx(0) = /275, 全 cx (24) 


由 式 (15) 及 式 (24) 可 获得 结构 地 震动 响应 的 
位 移 速度 及 层 间 位 移 及 其 变化 率 的 方差 。 
3.2 结构 随机 响应 的 功率 谱 分 析 

由 随机 振动 理论 ”"” ,结构 随机 响应 的 功率 谱 
密度 函数 与 其 协 方差 的 关系 为 


Si(w) = | Cir)eos(wr)dr (25) 


式 中 Sy(w) 为 结构 响应 XX 的 功率 谱 密度 函数 。 
把 式 (22) 代 入 式 (25) , 即 
2(n+2) 
Sx(w) = 25, > | ercos(wT)d7 (26) 
式 (26 ) 进行 积分 并 整理 , 即 
2(n+2) 
Pi 
Sx(w) 二 250 by Gx, ww +p? (27) 
观察 式 (27) ,结构 随机 地 震动 响应 的 功率 谱 表 
示 为 结构 复 振动 特征 值 与 频 域 变量 的 平方 和 的 倒数 
的 线性 组 合 , 为 结构 随机 地 震动 响应 0 阶 .1 阶 和 2 
阶 谱 矩 的 简明 解 葛 定 了 基础 。 
3.3 ”随机 响应 的 0~2 阶 谱 矩 分 析 
依据 平稳 响应 谱 类 的 定义 ”中 , 谱 矩 wx 表示 为 


ax = 2 | So)odw 


(28) 
式 中 71=0,1,2 分 别 表示 结构 随机 响应 的 0 阶 ,1 阶 
和 2 阶 谱 矩 。 
把 式 (27 ) 响应 功率 谱 代 入 式 (28 ) ,并 令 1=0， 
则 


2(n+2) 


. pi 
ts 三 下 5 cx | _ 
X,0 "之 瑟 4p 


对 式 (29) 积 分 , 则 结构 响应 的 0 阶 谱 矩 为 
ax0 = 2750 六 ou (30) 
由 平稳 随机 振动 理论 ,结构 地 震动 响应 的 0 阶 
谱 矩 2 阶 谱 算 与 结构 地 震动 响应 方差 的 关系 为 
(31a) 
(31b) 


dw (29) 


2 
Qx,0 = Ox 

2 
Qx,2 二 OCX 


式 中 X 表 示 X 对 时 间 的 一 阶 导 数 。 


动 响应 的 协 方差 表示 为 结构 振动 特征 值 指 函 数 的 线 
性 组 合 ,具有 表达 式 简 洁 明了 的 特点 。 

依据 随机 振动 理论 , 当 7 =0 时 ,结构 随机 地 震 
动 响应 的 协 方差 等 于 随机 响应 的 方差 , 即 


式 (31a) 表 示 建 筑 结构 的 随机 地 震动 响应 的 方 
差 与 式 (22) 表 示 的 方差 表达 式 完全 一 致 ,验证 了 本 
文 方法 的 正确 性 。 

将 式 (27) 表 示 的 功率 谱 代 入 式 (28) ,并 令 1=1 
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时 ， 倘 构 随机 地 震动 响应 的 1 阶 谱 矩 为 
3 idw (32 
a Ds 1 ee 人 (32) 


对 式 (32) 进行 积 4 


2(n+2) 


2(n+2) 2(n+2) 


= 250 limlne > Gxip; — 250 p> Gxipilnp; (33) 


由 文献 [10] 及 式 (28) 可 知 ， 2 阶 谱 矩 存在 , 则 1 
阶 谱 矩 肯定 存在 , 即 
> ip: =0 (34) 
由 于 实际 频 域 变 : 量 为 有 限量 ， 故 下 式 成 立 。 
2So limlno D2 p; =0 (35) 
F 则 结构 地 震动 响应 的 1 阶 席 逢 ax 可 表示 为 
0 a S28 > pilip? (36) 


(00) 及 式 (36) 可 和 结构 地 震动 响应 的 
谱 矩 均 有 解析 解 , 且 表 达 式 简洁 明了 。 


筑 结 构 体 系 动力 可 靠 度 


了 基 于 首次 超越 破坏 准则 的 建筑 吉 构 动力 可 靠 
ee 计算 是 线性 结构 可 靠 度 计算 的 最 重要 方法 。 
Vinarcke'” 提出 的 动力 可 靠 度 灾变 函数 的 理论 公 
5 询 具有 较 高 的 精度 ,得 到 了 广泛 应 用 。 

7) 
OY 2 1 -exp{ - 
ep| | 


TO x 


V2 (qx) Lr/ ox} 

s=0;3 (Lx/ox) ] 
(37) 

式 中 刀 为 响应 量 X 的 安全 界限 值 ;qx 为 谱 参 数 , 计 


算 公 式 为 
ax 1 
i (38) 
沪 Qx,0Qx,2 


则 构件 的 动力 可 靠 度 表示 为 
-Th(Zx)] (39) 
式 中 7 了 为 一 次 动 答 载 的 持续 时 间 。 

对 于 多 层 建 筑 结构 的 体系 可 靠 度 , 欧 进 萍 等 ” 
提出 了 “最 弱 失 效 模型 ”理论 ,其 以 层 间 变形 为 指标 
的 动力 可 靠 度 分 析 可 看 作 串 联 形式 的 结构 失效 模 
式 , 则 建筑 结构 体系 的 动力 失效 概率 P 为 


P,=1- ITP， (40) 
el 


20x) 1 -exp[ 


Py, exp[ 


5 算 例 应 用 


一 6 层 框架 结构 ,结构 层 的 力学 参数 :1 层 .2 层 
的 质量 均 为 520 t,3 层 ~6 层 质 量 为 450 t;1 层 ~2 
层 抗 侧 刚 度 为 1.8 x107 kN/m,3 层 ~6 层 抗 侧 刚 度 
为 1.0 x10 kN]m; 结 构 阻 尼 比 为 0.05;1 层 ~2 层 
层 高 为 4.8 m,3 层 ~6 层 层 高 为 4.0 m。 场地 地 震 
烈度 为 6 度 ,1 类 场 ,依据 文献 [6] ,地 震动 参数 取 值 
为 : w。 = 15.71 rad/s, &, = 0.72; S, = 11.87 x 
10 em /s’; w=0. 15w,,é€, =é,6 
5.1 本 文 解法 的 验证 

0 
度 功率 谱 和 结构 位 移 的 0 阶 、1 阶 和 2 阶 六 
S.1.1 ee 

由 式 (2a) 和 式 (2b) 得 

Wi = -2é,0,u, -Wis +2€ 0 %, + WIX, (41) 

根据 式 (6a) 式 (7) 及 式 (12) ,vs us xs sxs 可 

表示 为 


,分 别 对 地 面 加 速 
普 矩 开展 。 


t 


人 ua(t -7)dr (42a) 


i=1 0 


ktD 
未 
tD 


Xs = Uz = 


1 J Tr urn(t -7)dr (42b) 
0 


2(n+2) t 

i 本 ee -pir 

Us = Un+32 一 4 hs |e L urlt 7T) d7 
i 


© 


(42c) 


Ra ur(t T) d7 


(42d ) 
式 中 ,强度 系数 ;= wi, 本 算 例 的 强度 系数 见 
表 1。 
将 式 (和 2) 代 入 式 (41) 得 


pI 
> pe 
式 中 ,p; = - 26 /WA nr4 
wiAL, 3 
根据 式 (32) ,*。 的 单 边 功率 谱 密度 函数 为 


2(n+2)2(n+2) 
BB: pi 
S: (wm) = 45, 一 一 一 (44) 
by 和 2 pi +piw +p? 


ur(t — 7T)dr (43) 


骂 
| CrA2 十 2E WA 43 + 


表 1 , 模 态 强度 系数 表 


Tab.1 Modal strength coefficient of %, 
编号 复 特 征 值 特征 值 模 值 / 阻尼 比 强度 系数 
(rad.s ) 
1 1.697 二 1.635i 2.3565 0.72 1.7000 0.1070i 
2 0.2064.124i 4.1292 0.05 0 
3 0.572 二 11.431i 11.445 5 0.05 0 
4 11.311 干 10.902i1 15.7100 0.72 -13.011 于 2.453 9i 
5 -0.874 于 17.461i 17.4826 0.05 0 
6 -1.166+23.307i 23.3367 0.05 0 
7 —1.395+27.868i 27.9038 0.05 0 
8 -1.588 干 31.731i1 31.7712 0.05 0 


由 表 1 可 知 ,强度 系数 一 列 中 仅 模 态 1 和 模 态 4 
的 值 拓 零 , 且 模 态 1 及 模 态 4 的 振动 特征 值 的 模 值 
和 模 态 阻尼 比 与 C-P 谱 的 振动 特征 值 一 致 。 由 式 
(多 阿 知 , 半 仅 与 滤波 的 2 组 特征 值 有 关 , 而 与 建筑 
乡 梅 的 力学 参数 无 关 ,说 明了 本 文 方法 的 正确 性 。 
图 1 为 所 提 方 法 的 式 (44) 与 式 (4) 对 比 图 。2 
着 阴线 完全 吻合 ,验证 了 本 文 方法 的 正确 性 ,但 本 文 
罗湖 所 得 功率 谱 表 达 式 比较 简洁 。 


2.0r 
本 站- » 
» I 传统 方法 
wl5r A 一 本 文 方法 
> 
> 
SE, ) ~ 
S05 
= 泛 | 
| 
0 S100 L359 020 25 30 33 40 430 30 


wi(rad* s') 


图 1 5; (wm) 图 
Fig.1 Diagram of S; (&) 

5.1.2 结构 位 移 谱 短 

为 验证 所 提 计 算 谱 矩 方 法 的 正确 性 ,可 通过 式 
(30) . 式 (36) 对 结构 各 楼 层 位 移 的 0 ~2 阶 谱 矩 进 
行 计算 ,并 与 虚拟 激励 法 (PEM ) 进行 比 对 ,如 图 2 ~ 
图 4 所 示 。 由 于 PEM 法 在 计算 谱 和 矩 时 为 梯形 数值 
积分 法 , 需 给 定 频 域 变量 w 的 积分 区 间 和 积分 步 长 。 
本 文 的 积分 区 间 为 [0,1 000] ,积分 步 长 分 3 种 工 
况 : 人 四 工 况 1,Aw =0.8 rad/s;@ 工 沈 2,Aw =0.4 
rad/s;@) 工 况 3,Aw =0.25 rad/s。 

由 图 2 ~ 图 4 可 知 ,虚拟 激励 法 计算 谱 矩 随 着 积 
分 步 长 Aw 的 变 小 ,与 本 文 方法 越 来 越 接近 ,说 明了 
本 文 方法 的 正确 性 。 此 外 ,本 文 方法 的 CPU 耗 时 为 
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0.051 s, 虚 拟 激 励 法 工 况 1 的 耗 时 为 1.092 s, 工 况 2 
的 耗 时 为 2.081 s, 工 况 3 的 耗 时 为 5. 180 s。 本 文 方 
法 与 工 况 3 非常 接近 ,但 耗 时 仅 为 其 1/100 ,说 明了 
本 文 方法 的 高 效 性 。 
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Fig.2 Diagram of 0-order spectral moment of displacement 


of each layer of structure 
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图 3 结构 层 位移 1 阶 谱 矩 


Fig.3 Diagram of 1-order spectral moment of displacement 


of each layer of structure 


| 


图 4 结构 层 位 移 2 阶 谱 矩 


Fig.4 Diagram of 2-order spectral moments of 


displacement of each layer of structure 
5.2 结构 体系 动力 可 靠 度 

依据 抗震 规范 ,建筑 框架 结构 在 多 遇 地 震 下 的 
层 间 弹性 角 取 值 为 1/550 ,一 次 地 震动 持 时 为 15 s。 
利用 所 提 方 法 ,分 析 各 楼 层 层 间 位 移 的 0 ~ 2 阶 谱 
和 矩 . 谱 参数 及 动力 可 乱 度 , 见 表 2。 

由 表 2 可 知 ,各 楼 层 的 结构 动力 可 靠 度 均 较 小 ， 
按照 式 (44) 计算 的 结构 体系 失效 概率 为 1 ,说明 结 
构 整 体 失效 ,需要 增加 结构 抗 侧 刚度 或 者 设置 阻尼 
器 以 达到 减 震 的 效果 。 
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表 2 可 靠 度 分 析 表 
Tab.2 Reliability calculation of structure 
楼 层 0 阶 矩 / 1 阶 矩 / 2 阶 矩 / 构件 
a m2 (m2.s-1) (m2.s-2) 可 靠 度 
1 层 2.58x10-5 1.32x10-4 9.32x10-4 6.57x10-4 
2 层 ?2.18xl10- 1.03x10-4 6.39x10-4 8.93x10-3; 
3 层 5$.57x10- 2.52xl10-4 1.41x10-3 4.99x10-" 
4 层 4.10xl10- 2.01xl10-* 1.24x10-3 2.93x10-8 
5 层 2.53xl10- 1.56xl10-4 1.26x10-3 3.50x10-4 
6 层 9.16x10-56 7.41x10-5 8.21 x10-4 1.59x10-! 
人 、 
6 结 论 
下 A 
= 本文 针对 C-P 谱 激 励 下 建筑 结构 响应 和 动力 可 


HO 


于 旗 励 。 算 例 表 明 , 建 筑 结构 的 振动 特性 对 于 x, 
源 哆 响 ,说 明了 所 提 方法 仅 为 一 种 数学 方法 ,同样 适 
用 腔 其 他 复杂 的 非 白 噪声 平稳 地 震动 激励 模型 下 的 
结 攀 随机 地 震动 响应 分 析 。 


2) 所 提 方 法 所 获得 的 结构 地 震动 响应 方差 均 具 


解析 解 ,因此 本 文 方法 的 计算 效率 相 比 传统 方法 
:大 提高 ;同时 ,获得 了 结构 响应 的 0 ~ 2 阶 谱 答 
的 韦 析 表达 式 , 故 可 用 于 基于 Markov 分 布 的 可 靠 度 
分 析 ,进而 提高 了 可 靠 度 分 析 的 精度 。 

-三 3) 利 用 C-P 谱 的 滤波 方程 与 结构 运动 方程 组 成 
的 温 合 方程 为 非 经 典 阻尼 结构 ,采用 复 模 态 法 求解 ， 
因此 所 提 方 法 也 适用 于 复杂 的 且 具 有 线性 滤波 器 的 
随机 激励 下 各 类 线性 被 动 控制 结构 的 响应 计算 及 可 
靠 度 分 析 。 
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